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ABSTRACT 

We show that the factor ~2 discrepancy between spectroscopic and photo- 
metric age determinations of the Galactic bulge main-sequence turnoff can be 
naturally explained by positing an elevated helium enrichment for the bulge rel- 
ative to that assumed by standard isochrones. Helium-enhancement relative to 
standard isochrones is confirmed at the 2.3(T level. We obtain an upper bound on 
the helium enrichment for the metal-rich ([Fe/H]:^ +0.30) stars of AF ^ +0.11 
relative to canonical expectations, given the requirement that the spectroscopic 
and photometric ages be consistent and the limiting condition of instantaneous 
star formation. We discuss phenomenological evidence that the bulge may have 
had a chemical evolution that is distinct from the solar neighborhood in this 
manner, and we make several testable predictions. Should this emerging pic- 
ture of the bulge as helium-enhanced hold, it will require the development of 
new isochrones, new model atmospheres, and modified analysis and cosmological 
interpretation of the integrated hght of other bulges and elliptical galaxies. 

Subject headings: Galaxy: Bulge 



Introduction 



The mean age and age-spread of the Galactic bulge, teuige and Ateuige, are fundamental 
parameters of Galactic evolution, with hierarchical formation theories predicting an older 
bulge with more rapid formation timesca l e than models where bulge formation was linked t o 
early disk ( iKormendy fc Kennicuttll2004j : ICescutti fc Matteuccill201ll : llnoue fc Saitohll201ll ). 
However, age determinations of bulge stars are currently in a state of dissonance, with 
spectroscopic ages extending to much lower ages than photometric determinations. 

A seemingly secure picture had emerged from Hubble Space Telescope (HST) pho- 
tometry of the Galactic bulge main-sequence turnoff (MSTO). iKuijken fc RichI (|2002[ ) and 
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Zoccali et al.l ( 120031 ) both showed that the bulge MSTO is dim relative to the bulge hori- 
zontal branch (HB), with a brightness offset in distinguishable from that found in Galactic 
globular clusters (GCs), suggesting an old age. IClarkson et al.l ( 120111 ) refined this approach 
by searching for eclipsing binary blue-straggler stars near the MSTO. They found that what- 
ever sparse population of bright stars there was could be entirely explained by bulge blue 
stragglers, and concluded that no more than 3% of bulge stars are younger than 5 Gyr. 
Thus, the statement that the bulge stellar population is old appeared robust: bulge MSTO 
stars are dim. 



However, this picture has been challenged by lBensby et al.l ( 120101 . 1201 ll ). who used grav- 
itational microlensing events to obtain high-resolution spectra of bulge MSTO and subgiant 
(SGB) stars, which yielded detailed measurements of log (7, Teg, [Fe/H] and 12 additional 
chemical abundances. Their sample, microlensed source stars (only ~3% of microlensed 
sources are expected to be disk stars) with coordinates (|/| ^ 5,2 < |6| < 5) and radial ve- 
locity dispersion of ~100 km/s, is expected to be a predominantly bulge sample. Comparison 
of the measurements to isochrones led to a big surprise: a vast population of young stars. 
Seven of their 26 spectra (27%) yield ages, tinferrcd < 5 Gyr, an outcome with probability 
P < 10~^ if one assumes the 3% constraint. This method may seem complex at first glance, 
but upon close inspection it also appears robust. Consider the MSTO sta r MOA-2008-BLG - 
311S. At its high metallicity [Fe/H]= +0.36 ±0.07, standard isochrones dPotter et aPboOsI ) 
predict it should be no hotter than T^^ ^ 54507^ if one assumes an old age. However, its 
measured temperature is T^a = 5944 ±68i^, yielding the spectroscopic fit tinferrcd ~ 2.9 Gyr. 
Next, consider the metal-rich ([Fe/H]= +0.37 ± 0.05) SGB star MOA-2009-BLG-259S. At 
its temperature, Tcs = 4953 ± 93K, the expected surface gravity is logg ~ 3.90. However, 
the measured value is logg = 3.40 + 0.24, for tynfRrreH 3.0 Gyr. Therefore, the paradigm- 



challenging data reported by iBensby et al.l ( I2OIOI . 1201 ll ) pass the consistency test of being 
independently demonstrated in two distinct phases of stellar evolution. In summary: 
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Galactic GCs suggests teuig 



Brown et al. 


2010; 


Clarkson et al. 


2011 


). ' 


,ige ^ 12.8 Gyr ( 


Marin- Franch et al. 


2009) 



High-resolution spectroscopy of [Fe/H]>0 stars: hig h Tp^ for the MSTO and lo w logg 
for the SGB, indicating a young stellar population (IBensby et al.ll2010l . 1201 ll ). They 
find a weighted mean age of 5.9 Gyrs for the super-solar metallicity stars. 



In this Letter we suggest that the origin of this discrepancy is with the isochrones that 
are used to interpret the data and estimate the ages. Specifically, we question whether the 
adopted assumption of scaled-solar helium abundance is a valid approximation for the chemi- 
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cal evolution of the Galactic bulge at the metal-rich end. We demonstrate that the resulting 
impact on stellar parameter determination leads to a bias in inferred ages: photometric 
determinations are too old and spectroscopic determinations are too young. 
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Fig. 1. — Applying standard isochrones to helium-enhanced populations overestimates pho- 
tometric ages and underestimates spectroscopic ages. TOP: [Fe/H]=0, Y=0.274 (solar) 
isochrones with ages 6-13 Gyr in black (thicker for t=10 Gyr), and an [Fe/H]=0, Y=0.400 
(enhanced), t=10 Gyr isochrone plotted in red. On the spectroscopic plan, helium-enhanced 
populations mimic younger, standard populations with the effect minimized near the MSTO. 
BOTTOM: Same color scheme as above, but with tracks shown on the photometric {V — I, I) 
plane. The 10 Gyr helium-enhanced track has a much fainter MSTO, equivalent to a 
t~13.5 Gyr track of the same metallicity and standard helium. A large difference in he- 
lium. Ay =0.126, has been selected for demon strative purposes o f clarity. Isochrones are 
from the Dartmouth stellar evolution database (IDotter et al.ll2008l ). 
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2. The Degeneracy Between Age and Helium 



In general, isochrones assume the following prescription for initial helium: 



Y = Yp 



l-Yv 







1 + 



Ay 

AZ 



J0+ 



(2) 



where Y is the initial helium abundance of the star, Yp is the primordial helium abundance 
from big bang nucleosynthesis, Z is the metallicity mass fraction of the star, and (Ay/AZ)© 

i s the slope derived by fitting a line from the primordial values, ( Zp, Y-p) = (0, 0.249) 

(ISimha &: Steigmanll2008l ). to the solar values, (Zq, Yq) ^ (0.018, 0.272) ( Ivan Saders fc Pinsonneault 
201ll ). yielding {AY/AZ)q ~ 1.5. Some variation occurs as some isochrones use the pre- 
WMAP value Yp = 0.235. Equation ([1]) is manifestly problematic. There is no a priori 
reason that Y should be a single-parameter, deterministic function of Z; that such a func- 
tion should be a first-order polynomial; or that AY/AZ should be a universal constant. Why 
is the most abundant non-trivial element not generally incorporated as a variable input into 
isochrones? It is simply too difficult to measure: as the noble gas with two protons, it has 
the highest first ionization potential of any element, 24.6 eV. For the purposes of this work, 
we define a helium-enhanced stellar population as being a stellar population with values of 
Y exceeding those predicted by Equation ([T]). 

The effect of helium abundance on stellar evolution and thus inferred ages can be 
significant, due to the higher mean molecular weight and lower initial central hydrogen 
abundance accelerating the evolution of helium-enhanced stars. Hence, applying stan- 
dard isochrones to helium- enhanced populations leads one to overestimate photometric ages 
(IMarin-Franch et al.ll2010[ ). And further, as we will demonstrate, the opposite holds for 
spectroscopic determinations: ages will be underestimated. At fixed age and metallicity, 
MSTO stars with higher initial helium will be hotter and more compact, and thus bluer and 
dimmer. Their higher temperatures mimic younger ages on a log g-T^g diagram, but their 
smaller sizes and thus lower luminosity yield a dimmer MSTO, mimicking older ages on a 
CMD. On the SGB, the temperatures and bolometric luminosities become more similar, but 
not the surface gravities, causing spectroscopic ages to be underestimated. These effects are 
easily discerned in Figure [H 

The assumption that AY/AZ is a constant has been demonstrated to be a spectacular 
failure in GCs. For example, u Cen, NGC 2808, and 47 Tuc have a strong diversity of 
morphology in their CMDs, a fact that is we l l-explained by invoking helium-enriched second 



stellar generations f jPiotto et al.l |2005| . 120071 : ICarretta et al.l l2009l : Idi Criscienzo et al.l 12010 
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Nataf et al.ll2011by dl: iMilone et al.ll201ll : IValcarce fc Catelanll2011[ ). Though there is no ob- 
vious causal imphcation for the bulge, GCs do empirically demonstrate that AF/AZ can 
indeed vary. 



3. The Case for Elevated Helium Enrichment in the Bulge 



The helium abundance of the bulge had been investig ated in prior deca d es by compariri g 



the number o f red giant (RG) branch stars to HB stars. iTerndrupl (119881 ). iRenzinil ( ll994j ) 
and iMinnitil Jl995h respecti vely estimate d Y = 0.30 ± 0.05, Y = 0.30 - 0.35, and Y = 



0.28 ± 0.02. More recently, iNataf et al.l (l2011a| ]dl) found evidence for a helium-enhanced 
bulge using a different diagnostic: the red giant branch bump (RGBB). The bulge RGBB 
was shown to have anomalously low number counts (Aln(n) = —0.224 ± 0.080) and high 
luminosity (AJ = —0.104 ± 0.045) relative to expectations from Galactic GCs and stellar 
models, consistent with an enrichment AY ^ +0.06 at the median me tallicity. Both these 
methods may require incorporation of new evidence on RG evolution. iBrown et al.l (120111 ) 
demonstrated that 20-30% of low-mass white dwarfs are single, implying that a significant 
number of stars might skip the HB phase, suggesting the above work may have systematically 
underestimated the bulge helium abundance. 

The bulge should be expected t o have a higher AY/AZ than the Sun du e to it s 



elevated a-abundances ( Catelan 



2007h. iFulbright et al.l (120071). 



Alves-Brito et al. 



Bensbv et al.l (l2011f l. IConzalez et al.l (120111 1. Ijohnson et al.l (I2OIII ) and iRich et al.l (I2OIII ) 
each find a-abundance patterns elevated relative to the trends of the thin-disk, and possibly 
the thick-disk. The thin disk is the only stellar population for which the helium-metallicity re- 
lation has precise measurements (e.g., the Sun). Elevated [a/Fe] implies a lower relative con- 
tribution to chemical enrichment from type la SNe. Since these explosions have an effective 
Ay/AZ^ 0, reducing their relative contribution increases AF/AZ. Ad ditionally, a low con- 
tribu tion from type la SNe implies a rapid timescale for star formation (iHoriuchi &: Beacom 
20101). This in turn implies a higher relative contribution from the asymptotic giant branch 
(AGB) ejecta of higher mass stars, which are more helium-enriched than those of low and 
intermediate mass stars that take longer to reach the AGB. iRenzinil (120081 ) showed that the 
abundance of gas ejected from 3-10 Mq AGB stars is (Y) = 0.33, with the total yield and 
composition being a function of the metallicity and the slope of the initial mass function. 
Further, stellar models predict that the winds of fast rotatin g massive stars (WFR MS) are 
extremely helium-enriched, with values reaching (Y) = 0.50 ( iDecressin et al.l 120071 ). Thus, 
the a-enhancement of the bulge implies a higher effective value of AY/AZ due to two dis- 
tinct causes. As the duration of star formation increases, lower-mass AGB ejecta, which 
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are less He-enriched, will contribute relatively more to the chemical evolution budget than 
those of higher mass AGB+WFRMS ejecta. Second, late-time Type la SNe, explosio n s that 
contribute a lot of metals but no helium, will begin to go off. Figure 1 of ICatelanI (120071 ) 
implies Yf»0.31 at [Fe/H]=0 for a population with an a-enhancement similar to the bulge . 
This is midway between standard isochrones values and the estimate of iNataf et al.l (j2011dl ). 
We note that there has been some skepticism regarding WFRMS as the source of GC helium 
because they enrich metals as well as helium, and because GCs have a weak gravitational 
potential well that may not hol d on to these ejecta . Neither of those concerns apply to 
the bulge. The abundances from iBensby et al.l (120111 ) suggest that the sodium-oxygen anti- 
correlation seen in GCs and thought to trace high-mass AGE and/or WFRMS ejecta, is 
present in metal-rich bulge stars and trends with [Fe/H]. 



We refer to a new class of models (IChung et al.l 1201 ll : iBekkil 120121 ) motivated by the 
need to explain the UV-upturn of old bul ges and ellipti c als, w hich are bes t expla ined by the 
presence of blue horizontal branch stars. IChung et al.l (120111 ) and iBekkil (120121 ) both posit 
imperfect mixing between AGB and SNe ejecta, thereby having stars form directly from the 
gas of helium-rich AGB ejecta without dilution from SNe ejecta. If this process occured in 
other spheroids, it likely also occured in the bulge, and vice versa. These models also make 
the testable prediction that there should be a large helium spread at fixed metallicity. 

There is another eff ect that may be a t play : helium sedimentation. Simulations of hot 
gas in galaxy clusters by IChuzhoy &: Loebl (12004 ) predict that helium preferentially accretes 
onto cluster cores, as: 



Ay 



t 



10 Gyr/ VO.l 



gas 



T \3/2 

10 Kev) 



(3) 



They argue that values of AY ^0.12 can be reached by cD galaxies, which could explain 
the UV upturn of these galaxies by creating more extreme blue HB stars. While the effect 
in the bulge should be much smaller, it still may contribute to helium enhancement. 
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Evidence of Enhanced Helium in the Analysis of Spectroscopic Data and 

Predictions Thereof 



We discuss two systematic trends predicted by stellar models. First, we predict that 
the age offset induced by applying standard spectroscopic {logg — T^s) isochrones should be 
minimized near the MSTO and maximized on the SGB: the highest spectroscopic ages should 
be for stars nearest to the MSTO. Second, we predict that the error in the inferred absolute 
magnitude (M/ true — M/ inferred) should on average be positive, and correlate p ositively with 
the sp ectroscopically inferred mass. The predicted trends as well as the data of iBensby et al. 
( 1201 ll ) are shown in Figures |2] and El Both predicted trends are consequences of the same 
underlying phenomenon; the size of the age offset being a function of evolutionary state. 
The theoretical predictions shown in Figure [2] are summarized in Table [TJ 



We use isochrones from the Dartmouth stellar evolution database (jPotter et al.ll2008l ). 
The 8 input stellar parameters are {t =10,13 Gyr, [Fe/H]=0, [a/Fe]=0,+0.4, 0.33,0.40), 
where ( [Fe/H]=0, [Qi/Fe]= +0.40) is approximately equivalent to ([Fe/H]= +0.30, [a/Fe] = 
+0.00) (ISalaris et al.lll993l ). The upper four panels in each figure show the effect of different 
offsets in helium: Ay =0.028, 0.056, 0.098 and 0.126, for three assumed intrinsic ages, 6 Gyr 
(magenta), 10 Gyr (blue) and 13 Gyr (red). We evaluate these populations using standard 
isochrones {Y = 0.245 + 1.6Z) from the same database. When analyzing predicted stellar 
observables of one age and composition with isochrones of a different age and composition, 
we use the metric: 

Alog^ 



AS' 





2 


[92 kJ 


+ 



0.195 dexJ 



(4) 



where the normalizations are the average errors reported by lBensby et al.l (120111 ). Inspection 
of Figures [2] and [3] reveals that the predicted effects are monotonic with AF at fixed log g or 

(M/Mq) 

Inferred • 



- 8 - 




Fig. 2. — Applying standard isochrones to helium-enhanced populations leads to an offset 
between the inferred and true age that is monotonic in Ay and is minimized near the MSTO. 
Top four panels show the inferred age for stars 6 Gyr old (magenta), 10 Gyr old (blue) and 
13 Gyr (red) by standard isochrones as a function of logg, for different values of helium- 
enhancement. The true ages are shown by the horizontal lines to guide the ey e. Bottom 
panels shows the reported ages as a function of log g from iBensby et al.l (120111 ) , separated 
into subsolar and supersolar samples. The metal-poor stars have a flat log g-age distribution, 
whereas the metal-rich stars show older ages near the MSTO. 
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Fig. 3. — Top four panels show the predicted error in the absolute magnitude inferred for 
stars 6 Gyr old (green), 10 Gyr old (blue) and 13 Gyr (red) by standard isochrones of the 
correct metallicity as a function of the inferred stellar mass, for different values of helium- 
enhancement. Bottom panels show the difference between absolute magnitude inferred from 
the microlensing lightcurves and th e inferred absolute m agnitude as a function of inferred 
spectroscopic mass for the sample of iBensby et al.l (120111 ). The best- fit line to the metal-rich 
stars has both a greater intercept and a greater slope. 
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5. Preliminary constraints of the Age and Helium Abundance of the Bulge 

We first estimate Ay given the assumption that isochrones with the "correct" helium 
values would yield ages on the MSTO and SGB that are independent of logg. We test for 
this by rescaling the ages (and their errors) using the values of tinferred/^true from Table [H 
and we compute x^: 

We obtain = 30.3, 28.2, 27.2, 24.8 and 25.6 for AY = 0.0, 0.028, 0.056, 0.098, and 0.126. 
Thus is minimized for AY = 0.098, yielding Ax"^ = 5.4 - a 2.3a detection of helium- 
enhancement. The inferred enhancement is therefore AY = 0.098 ±0.043 at [Fe/H] ^ +0.3. 
In contrast, applying the same procedure to metal-poor stars increases x^ . 

We now estimate, using a different method, an upper-bound on AY. We first impose the 
requirement that the photometric and spectroscopic ages agree. From Figure [T], we see that 
a 10 Gyr isochrone with AY = +0.126 has an MSTO as faint as that of a standard ~13.5 
Gyr isochrone, for a 35% offset, thus the photometric age error can be linearly approximated 
by: 

^fl + 2.8AFy (6) 

\ ^True ^ Phot \ / 

For the spectroscopic age inference, the age offset is a function of evolutionary state, so we 
estimate by taking the mean value of tinfcrred/i^Tme in the range 3.44 < logg < 4.60: 



/ ^inferred \ 

\ ^True ' Spec \ / 



2.9Ay), (7) 



and from the literature (IBensby et al.ll201ll : IClarkson et al.ll201ll ). we have that: 



''^Photometric < 2 Q " ^ ~^ 2.8AY 



^Spectroscopic 1 2.QAY 

Combining these equations yields a helium enrichment relative to isochrones of AY < 0.11 
for the metal-rich stars, yielding a corresponding mean age of t^uige ~ 10 Gyr, and a value 
of (Ay/AZ)Buige ^ 5.0 if one assumes a linear form for Y{Z). As the canonical helium 
abundance for stars at [Fe/H]Ri +0.30 is y ~ 0.31, this sets an upper bound of F ~ 0.42. 
This upper bound on AY is derived by equating the ages of the metal-rich bulge stars with 
the mean age of the bulge: it assumes the limiting condition of instantaneous s tar formation. 



Any a ge spread will decrease the amount of extra helium required. In contrast, iBensby et al. 



( 120111 ) derive an age difference between the metal-poor and metal-rich stars of ~6 Gyr, under 



the assumption of standard helium enrichment. 
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Discussion 



We have shown in this Letter that the discrepancy in Galactic bulge turnoff age estimates 
can be well-explained by the hypothesis that the chemical evolution of the Galactic bulge is 
helium-enhanced relative to that assumed by standard isochrones. We first found evidence 
for helium-enhancement at the 2.3a level, suggesting AY = -1-0.098 for the metal-rich stars. 
We also derived an upper bound on the helium enrichment of the bulge, AY < 0.11 for 
the metal-rich ([Fe/H];=s -1-0.30) stars. We deliver two testable predictions, that the absolute 
magnitudes will be systematically fainter than the best-fit values from standard isochrones, 
and that inferred ages will be highest near the MSTO. Both diagnostics suggest helium- 
enrichment for the metal-rich stars, with no such evidence for the metal-poor stars. 

There are multiple roads forward. On the theory side, chemical evolution models will 
be required to ascertain what conditions are required in the primeval Milky Way for a 
high value of (AF/AZ) Bulge- A broader range of isochrones than currently available will 
be required to interpret the available photometric and spectroscopic data. Finally, it is 
possible that the stellar atmosphere models used by spectroscopists fail in the case of at- 
mosp heres that are extre mely helium-enriched. The [Fe/H] distributi on of MSTO+SGB 



Hill et al. 



2011 



stars (jBensbv et al.ll201l|) is inconsis tent with that derived from giants (jZoccali et al.l 12008 



Johnson et al.l 120111 ). as it has a trough near [Fe/H]=0, where the giant 



metallicities show a peak. Could this be the reason? Observationally, an uncertainty in 
the value of (AF/AZ)Buige should be incorporated in analyses of bulge stars. The age dis- 
crepancy discussed in this paper constitutes the fourth, independent line of evidence for 
a helium-enriched bulge, following previous investigations of star counts on the upper RG 
branch, on the RGBB, and chemical evolution arguments based on the high values of [a/Fe]. 

The helium abundances suggested here are d ifficult to attain with chemical evolution 



models, as AYAZ < 3.0 for general assumptions (IChiappini et al.ll2002l : ICarigi fc Peimbert 



20081 ) . Should the inferred helium-enhancement for the bulge hold up under further scrutiny 
it may necessitate rethinking the underlying assumptions of chemical evolution models. 
Additionally, even in the disk, there may be surprises with respect to helium-enrichment: 
Gennaro et al.l ( 120101 ) find AY/AZ = 5.3 ± 1.4 in a sample of local K-dwarfs. 



The preliminary estimates of this Letter result from the assumption of instantaneous 
star formation. How much age spread is reasonable depends on assumptions of the initial 
condi tions of the Galaxy. The as sumption of rapid gravitational infall yields Atsuige ~ 0.5 
Gyr (ICescutti fc Matteuccil l201ll ). whereas the "clump-origin bulge" model, in which the 
bulge is formed by stellar clumps forming in a gas-rich disk and accr eting to its center 
by dynamical friction, yields A^Buige ~ 2.0 Gyr ( llnoue fc Saitohl l201ll ). However, recent 
investigations of star counts and radial velocities find that bulge kinematics are consistent 
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with a pure N-body bar that evolved from secular instabihties, which would imply that bulg e 



stars are just disk stars on bar orbits ( Shen et aLKOlulMartinez-Valpuesta fc Gerhardll201ll ). 



This would suggest a duration of star formation as extended as that of the inner disk. 

The results of these endeavors may have significant cosmological implications. Popula- 
tion synthesis studies of other galaxies assume a solar value of AF/AZ. If this parameter is 
demonstrated to be non-universal for field stars, it will affect age and mass determinations 
of field galaxies. 
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Table 1: Predicted values of tinferred/^TVue as a function of log(5f) and AF. The first four 
columns are for the ^True = 10 Gyr case, and the last four columns are for the tivue ~ 13 
Gyr case. 
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